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Employing a novel type of non-commutative product in the Dirac-Kahler twisted 
superspace on a lattice, we formulate a field theoretically rigid framework of extended 
supersymmetry on a lattice. As a first example of this treatment, we calculate one-loop 
■ (in some cases any loops) quantum corrections for a twisted Wess-Zumino model with 

N x N hermitian matrix superfields on a two dimensional lattice. The calculations 
are entirely given in a lattice superfield framework. We report that the mass and 
the coupling constant are exactly protected from the radiative corrections at non-zero 
lattice spacing as far as the planar diagrams are concerned, which implies the realization 
of exact lattice supersymmetry w.r.t. all the supercharges in the large- N limit. 



1 Introduction 



in 

There have been a variety of studies addressing the subject of supersymmetry (SUSY) 
on a lattice. The literal purpose of lattice SUSY is to provide a well-defined constructive 
framework of SUSY formulations with which one may extract the fully non-perturbative 



in 

information either analytically or numerically. A successful formulation of lattice SUSY 
may be expected to provide a deeper understanding of the relationship between bosons and 
fermions in a regularized framework. In recent years, it has been recognized that the so- 
called twisted SUSY formulations are playing particularly important roles in this challenging 
subject [U El [3, HI El [6] . It was also recently pointed out that amongst the other formulations 
of lattice SUSY [7] , the deconstruction formulation [8] is quite closely related to the twisted 
SUSY formulation on the lattice [9]. In spite of these developments, however, realizing all 
the SUSY generators in a SUSY algebra being considered, say, D = 4Af=lorD = 2Af=2 
or whatever @, exactly on a lattice, has been considered as a very non-trivial task. This is 
essentially because in general only the finite subgroup of the entire translational symmetry 
group in the continuum spacetime can be survived at non-zero lattice spacings. 

In the series of papers [TJ [21 [3], we have proposed an extended SUSY formulation on 
a lattice which keeps all the SUSY transformations exact at least at classical level. The 
formulation is based on the two crucial theoretical ingredients. The one is the notion of 
twisted SUSY. This is essentially originated from the intrinsic relation between twisted 
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fermions and staggered fermions or more appropriately Dirac-Kahler fermions [10J. It has 
been pointed out that the fermionic components in the lattice SUSY multiplet are realized 
as a staggered fermion, and accordingly the spectrum doubling turns to have a physically 
relevant meaning as an extended SUSY degrees of freedom pp. The lattice fermions in this 
framework have their own geometrical meaning. They have one-to-one correspondence with 
all the simplices available on a lattice being considered, namely, 0, 1 and 2- form on a two 
dimensional lattice or 0, 1, 2, 3 and 4- form on a four dimensional lattice, which is a natural 
consequence of the Dirac-Kahler fermion realization on a lattice [IT] . The other crucial 
ingredient, which is also tightly related to the Dirac-Kahler structure, is the "mild" non- 
commutativity associated with the lattice supercharges. It has been long recognized that 
the breakdown of Leibniz rule is one of the major obstacles in realizing exact SUSY on a 
lattice. The notion of "mild" non-commutativity has been introduced in order to overcome 
this difficulty The non- commutative nature of lattice supercharges can be trace back to the 
intrinsic non- commutative property of difference operators. In [1] , we have shown that if the 
non-commutativity matching conditions, which we refer to as the Leibniz rule conditions on 
a lattice, are satisfied, then all the supercharges can be exactly realized on the lattice at least 
at classical level. This notion of "mild" non-commutativity has been further applied to the 
formulations of extended supersymmetric gauge theories by means of link supercharges and 
link component fields [21 E] and recently also applied to the formulation of Chern-Simons 
theory on a three dimensional lattice [12J. 

In spite of these developments, on the other hand, we have not addressed the quantum 
aspects of these formulations so far since it has not been fully established how to deal with 
the above "mild" non-commutativity in a quantum field theoretically rigid manner. In this 
paper, we attack this problem and propose a field theoretically well-defined treatment of 
the "mild" non-commutative SUSY formulation on a lattice. In particular, we shall employ 
a novel type of non-commutative product in the Dirac-Kahler twisted superspace which 
honestly accommodates the "mild" non-commutativity introduced in Ref. [Tj. As a first 
example of this framework, we shall calculate one-loop (in some cases any loops) quantum 
corrections for a twisted Wess-Zumino model with N x N hermitian matrix superfields on a 
two dimensional lattice. The calculations are entirely given in the lattice superfield method. 
We report that the mass and the coupling constant are exactly protected from the quantum 
corrections as far as the planar diagrams are concerned, while the SUSY spoiling non-planar 
contributions are suppressed at least by 0(a 2 /N). These features imply the realization of 
exact lattice supersymmetry w.r.t. all the supercharges in the 't Hooft large- iV limit. 

This paper is organized as follows. In Sec. 2, after introducing a novel type of non- 
commutative product which honestly represents the "mild" non-commutativity which we 
have introduced in [TJ, we (re-)formulate the lattice counterpart of extended SUSY algebra 
and see how to construct an extended SUSY invariant action on a lattice. In Sec. 3, as 
a simple and explicit example, we introduce the Dirac-Kahler twisted J\f — D — 2 SUSY 
algebra and explicitly see how the lattice SUSY algebra is realized in accordance with the 
lattice Leibniz rule conditions. We then construct a M = D = 2 twisted Wess-Zumino 
model on a lattice and proceed to derive the superfield propagators and the three-point 
vertex functions, taking a full advantage of the lattice superfield method. In Sec. 4, we 
extend the chiral and anti-chiral superfields introduced in Sec. 3 to N x N hermitian matrix 
valued superfields in order to set up the Af = D = 2 twisted Wess-Zumino model with 
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global U(N) symmetry. We perturbatively investigate the possible radiative corrections to 
this model by explicitly calculating the one-loop (in some cases any loops) contributions to 
the kinetic term, the mass term and the interaction term of the lattice action. We then 
find out that the quantum corrections strictly protect the lattice SUSY as far as the planar 
diagrams are concerned. In particular, we shall see that apart from the overall wave function 
renormalization the mass and the coupling constant are strictly protected (at least to this 
order) from the radiative corrections in the planar diagrams. These features imply the full 
realization of exact SUSY at non-zero lattice spacing in the 't Hooft large- N limit. Sec. 5 
summarizes the formulation with some discussions. 

2 Introducing a Non- Commutative Product 

In this section, we first introduce a novel non-commutative product in the Dirac-Kahler 
twisted superspace which honestly represent the "mild" non-commutativity introduced in 
Ref. pp and present how to formulate a lattice counterpart of SUSY algebra. The non- 
commutativity introduced in Ref (Tj in order to keep the lattice Leibniz rule is the one 
between the bosonic (discrete) coordinate x^ and Grassmann coordinate 9 A in the twisted 
superspace (x^,9 A ), 

e A f(x) = {-l)Wf(x-2a A )0 A , (2.1) 

or equivalently, 

[x^Ba] = 2(o4) M #4, (2.2) 

where the symbol |/| takes value of +1 or —1 for the bosonic or fermionic function /, respec- 
tively. 9 a is representing any fermionic coordinate in the Dirac-Kahler twisted superspace, 
namely, 9 A = (9, 9^ 9) with p = 1, 2 in N = D = 2 pQ H, 9 A = (9, 9^9^,9) with fx = 1 ~ 3 
inA/" = 4D = 3 [3] and 9 A = (9, 9^ 0^, 9^ 9) with v = 1 ~ 4 in M = D = 4 [2]. The 
vector a a in the r.h.s. of (12. 2p governing the non-commutativity between x^ and 9 a is to be 
constrained at a later stage such that the resulting SUSY algebra be consistent with the non- 
commutativity matching conditions which we refer to as the lattice Leibniz rule conditions. 
In this framework, superfields may generically be expanded in terms of the non-commutative 
9 as such as 

9 a) = 4>{x) + 9 a <Pa{x + a A ) + 9 a 9 b <Pab{x + a A + a B ) H , (2.3) 

= <p(x) + (-l)^<p A (x-a A )9 A + cf> AB (x-a A -a B )9 A 9 B + --- , (2.4) 

where the component fields are chosen to be located in a symmetric manner. The symbol |0| 
again takes the value +1(— 1) for the bosonic (fermionic) component field <p. All the superfield 
operators such as supercharges and supercovariant derivatives are accordingly represented 
in terms of the non-commutative #a's on the lattice. The highlight of the formulation based 
on the relation (12.21) is that the resulting lattice action can be manifestly invariant under all 
the SUSY variations at least at tree level if the lattice Leibniz rule conditions are satisfied 
[I]. On the other hand, it is rather non-trivial how to deal with the non-commutativity 
(12.21) within the framework of quantum field theory. One possible way to provide a field 
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theoretically more rigid framework is to formulate the bosonic and fermionic coordinates 
entirely in terms of matrix representations and deal with the (super)fields as the functions of 
these (super) coordinate matrices. A formulation presented in Ref. [6J may be placed along 
this direction. 

A somewhat alternative way to treat the non-commutativity (12 .2p . which is our starting 
point of this paper, is to alter the notion of product and introduce a non-commutative 
product in such a way to satisfy the relation, 

[x^9 A ]* = x fl *9 A -9 A *x fl = 2(a A ) fl 9 A . (2.5) 

This type of non-commutativity treatment may be regarded as a certain analogue of the 
Moyal star product whose quantum implications have been extensively investigated (e.g. 
[13] for a review). The principle we shall take is that all the non-commutativity should 
be originated from the nature of the product denoted as * in (12.51) . while the x^ and 9 A 
should be treated as the ordinary (discrete) bosonic coordinates and fermionic coordinates, 
respectively. In this framework, superfields can be expanded in terms of ordinary fermionic 
coordinates and ordinary functions, |f| 

$(x,0 A ) = (f>(x) + 9 A( f> A (x) + 9 A 9 B( f> AB (x) + ■ ■ ■ , (2.6) 

and thus may be served as path integration variables such as in the standard SUSY treat- 
ment in the continuum spacetime [T4T [TBI [TBI [T7] . A star product of superfields, which is 
intrinsically non-commutative, is defined in terms of * subject to the relation (12.51) . 

A ) * <t> 2 (x, 9 B ) ^ $ 2 (x, 9 B ) * $ a (ar, 9 A ). (2.7) 

An explicit representation of the star product satisfying the relation (12.51) may be given by 

$i(z, 6 A ) * $ 2 (x, 9 B ) = ii(T L x T R $i(x, 6 A ) ® $ 2 (x, 9 B )), (2.8) 

with the elements T L and J-'r, 

T L = exp{V y^MpQAj^- ® dp], T R = &cp{y2d p ®y2(a A ) p 6 A -—}, (2.9) 

p A A p A d ° A 



where in (12. 8p the symbol \x denotes the multiplication map, //(/ <S> g) = fg. The operator 

a 



Oat&- m the r.h.s.'s of (12.91) should be understood as an adjoint operation 



^ w A adi " A W A = **k' (2 ' 10) 



together with the anti-commutation relation between and 9 B , 



{-^-,6b} = S ab . (2.11) 



3 We shall take a convention to locate all the component fields on a same site. Note that superfields can 
be re-expressed in terms of the star product. Using (|2. 14|) for instance, one has <&(x, A ) = 4>(x)+6a*4>a{x + 
o-a) + 0a * Ob * 4>ab{x + a A + a B ) H = <f>(x) + (-l)^- 41 ^^ - a A ) * A + 4>ab(x — a A — a B ) *0a*Qb-\ • 
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One may easily verify that the elements T L 1 and Tr independently give rise to the following 
relations, 

y,{TZ x 9 A ®f[x)) = 9 A f(x-a A ), M^I 1 f(x) ® 9 A ) = f(x)9 A , (2.12) 
^ R 9 A ®f{x)) = 9 A f(x), ti(F R f(x)®9 A ) = f(x + a A )9 A . (2.13) 

Accordingly, the star product defined in (12.81) would give, 

9 A * fix) = 9 A f(x-a A ), f(x)*9 A = f(x + a A )9 A . (2.14) 

It is now obvious that the relation ( 12. 5p is satisfied just by taking f(x) = x. It is important 
to notice that the star product ( 12. 8ft satisfies the associativity, 



{$i(x,e A )*$2(x,e B ))*Mx,9c) = $i(x,9 a )*(Mx,0b)*Mx,8c))- (2.15) 

d 
08a 



Also note that the star product of the super derivative operators gives rise to the opposite 



shifts compared to (12.14[) . 



d d d d 

f( x ) = T^rf(x + a A ), f( x )*^7T = f(x-a A )—-, (2.16) 



d9 A JV ' d9 A y n JX ' 09 A Jy J d9,_ 

with which one can verify that the anti-commutation relation ( 12. lip holds also in the star 
product language, 

r\ r\ r\ 

W 6b} * = w a * 9b * f{x) + dB *w A * f{x) = SABf ^- ( 2 - 17 ) 

Although the combination of the elements (!F£ ) t (J 7 R) 2 ~ t with an arbitrary parameter t 
would always induce a non-commutative product which satisfies the relation ( 12.51) . in this 
paper we shall only consider the symmetric definition given in ( 12.81) . 

Having introduced a non- commutative product between the bosonic and fermionic coor- 
dinate of the superspace, x^ and 9 A , we also require the same amount of non-commutativity 
between x M and Grassmann parameters £ A , 

U*f(x) = Uf(x-a A ), f(x)*U = f(x + a A )U, ( 2 -!8) 

such that the operation of ^Agf - would properly induce a derivative w.r.t. 9 A also in the 
language of the star product, 



d 

U-^-A{x) + 9 B (j)B{x) + 



d 

6»oT*Wi) + Mb(i) + - 

OU A 

d 

-(<P(x) + O b <I>b{x) + -~) *U- 



09 A 

= SabUMx) + " • • • (2.19) 

This extension can easily be done by just replacing the #Agfj as £Ag|j + 9 A -^ in the 
expressions ( 12.91) . 
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Based on the above setup, we now proceed to define a lattice counterpart of SUSY 
algebra. Since in the continuum spacetime, SUSY algebra can be regarded as a fermionic 
decomposition of differential operators, one may naturally expect that a lattice counterpart of 
SUSY algebra can be formally expressed as a fermionic decomposition of difference operators, 

{Qa, Qb] = /Vm, (2-20) 

where f^ B denotes a numerical coefficient, while denotes a "formal" difference operator 
to be specified in the following. A SUSY transformation on the lattice is defined by a * com- 
mutator between a supercharge Qa associated with a corresponding Grassmann parameter 
U and a superfield <3>(x, 6 C ), 

5 A <S>(x,6c) = [UQa,$(x,0 c )U 

= UQa * $0, C ) - $0, 9 C ) * UQa {A: no sum). (2.21) 

A commutator of 5a and 5 B would accordingly give, 

[5 A , 6 b ]*{x, e c ) = [UQa, [UQb, ®(x, 6 c )U* ~ [UQb, [UQa, 6 c )U* 

= -[UU{Qa,Qb},Hx,6 c )1 

= -UUfAB[ d ^ ( y x ~ a A-dB,Oc)-^{.x + aA + a B ,e c )d ll ], (2.22) 

where in the last line we substitute the relation (12.201) . One can notice here that the SUSY 
variation defined in terms of the star product (12.211) automatically induces a difference op- 
eration in the last line, provided either of the following conditions holds, 

a A + a B = +n fl and d M = — ^, (2.23) 

aA + a,B = —Tin and = +-, (2.24) 

where n M is denoting a unit vector whose component is defined by {n^) v = 5^ u . One may 
verify that the last line of the expression (12.221) has the desired continuum limit if the 
condition either f!2.23j) or (12.241) is satisfied, 

-UUf^l[Hx + n„9c)-Hx-n„e c )} - -&£b/Wm*(Mc)- (2.25) 

The superspace expressions for the supercharges satisfying (I2.20|) are readily given by, for 
instance, 

Qa = -^- A +\fA B dBd„ Q B = ^ + \fAB0Ad,. (2.26) 

Remark here that the difference operation in this framework is generated by the non- 
commutative nature of the star product and thus the lattice SUSY algebra given in (12.201) is 
actually describing a "reduced" algebra in the sense that the differential operator <9 M in the 
continuum SUSY algebra is replaced by a constant value d^ = ±| in its lattice counter- 
part. The conditions (I2.23|) and (I2.24[) coincide with the Leibniz rule conditions introduced 
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in the previous works [H [2l E] where it was shown that Dirac-Kahler twisted SUSY algebra 
of Af = D = 2, Af = 4D = 3 and J\f = D = 4 can satisfy these conditions. Although the 
above star product formulation is formally available in any of these cases, in this paper we 
mainly concentrate on the simplest but non-trivial case, M — D — 2. It is also important to 
remark that a SUSY variation for a star product of superfields $i * $2 satisfies the following 
Leibniz rule thanks to the associativity (e.g. (12.151) ). 

6 A ($i * $2) = [UQai $1 * $2]* 

= UQA * $1 * $ 2 - $ x * $ 2 * £ A Q A 

= (5 A ®1 + $1 * UQA) * $ 2 - $x * (UQA * $2 - Sa®2) 

= 5 A &i * $ 2 + $1 * 6 A ®2, (2.27) 
while a SUSY variation for an ordinary product $i$2 does not, 

<5 A ($!$ 2 ) 7^ $a®i®2 + $1«A$!, (2-28) 

since in general we have, (A*B)C 7^ A* (BC). From these properties, one can immediately 
see that any lattice action which is to be invariant under the SUSY variation ( 12.211) should 
be expressed in terms of the star product. 

3 N = D = 2 Dirac-Kahler Twisted SUSY Algebra & 
Wess-Zumino Model 

In this section, as a simple and explicit example, we introduce the M = D = 2 twisted 
SUSY algebra and see how the lattice counterpart of the algebra can be realized in accordance 
with the lattice Leibniz rule conditions. We introduce a chiral and anti-chiral superfields to 
construct a J\f = D = 2 twisted Wess-Zumino model on a lattice. The superfield propaga- 
tors and the three-point vertex functions shall be explicitly derived by means of the lattice 
superfield method. 

3.1 M = D = 2 twisted SUSY algebra on a lattice 

Let us first introduce the M = D = 2 SUSY algebra in the continuum spacetime which 
is given by 

{Qd,Qjf>} = ZiSijdidapd^, (3.1) 

% % 
[J,Qad] = -A.lh)apQpii [R,Qai\ = ^(lb)ijQaj, (3.2) 

[JM = ie^d,,, [R,d„] = [d^Q ai \ = [d^d u ] = [J,R] = 0. (3.3) 

where the indices a, (3 = 1 or 2 and i,j = 1 or 2 denote spin and M = 2 internal indices, 
respectively. The gamma matrices in two dimensions can be represented by Pauli matrices, 

7i = °3, 72 = en, 75 = 7i72 = i&2- (3.4) 
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The conjugate supercharge Q = can be taken as Majorana, subject to the condition 
Qia = (Q T )ia = Qai- J an d R denotes the generators for the SO (2) Lorentz and 5*0(2) 
internal rotations, respectively. The twisted form of the algebra (13.11) is given in terms of 
the Dirac-Kahler expansion of the supercharge, 

Q ai = (1Q + 7mQm + 75<5)m, (3.5) 

as 

{Q,Qv} = id v., {Q,Qii} = -itfiudv, {others} = 0. (3.6) 

The twisted supercharges (Q, Q M , Q) transform as a scalar, a vector and a (pseudo-) scalar 
under the twisted Lorentz generator J' = J + R, respectively, in the continuum spacetime, 

[J',Q] = 0, [J',Q M ] = ie^Q„, [J',Q] = 0. (3.7) 

It was first pointed out in Ref. [T] that the algebra which can be realized on the lattice 
is not the form of (13. ip but the form of (I3.6P since only the Dirac-Kahler twisted form of 
SUSY algebra has a chance to satisfy the Leibniz rule condition. In the present context of 
the formulation, the lattice counterpart of the algebra ( 13.61) can be expressed as, 

{Q,Qii} = id+, {Q,Qn} = -ie^ u d~, {others} = 0, (3.8) 

with the corresponding Leibniz rule conditions, 

a + <V = +n^, a + a^ = -\e^\n u . (3.9) 

The "formal" difference operators and d~ in the r.h.s. of ( 13.61) take the value of — ~ and 
+|, respectively, as is required by the conditions (I2.23P and (12.241) . The conditions (I3.9p are 
satisfied by the following generic solutions with the one vector arbitrariness, 

a = (arbitrary), a^ = — a, a = — ri\ — n 2 + a. (3.10) 

Note that the sum of all the shift parameters vanishes regardless of the one vector arbitrari- 
ness, 

a + a + a\ + a 2 = 0. (3-H) 

In this paper, mainly because of the advantages in dealing with the superfield formulation, 
we shall take the symmetric choice which is depicted in Fig. [TJ 

a = (+ 2' + 2 ) ' 01 = (+ 2'~2 ) ' ° 2 = ( ~2' + 2 ) ' a = ( ~2'~2 ) - ( } 

Superspace expressions for the supercharges Qa = (Q, Q^, Q) which satisfy (13.81) are given 
by, 

Q = qo + Z 2 9 ^> Q » = W +l 2 dd »~ l 2 h » vd ~ v > (3 ' 13) 

— fj, — 

Q = is- far^ • (3-14) 





■ x + n-2 




a 


~ ni / 




a 


CLl 




■x - n 2 



Fig. 1: Symmetric choice of a a 



We also introduce the following supercovariant derivatives Da = (D, D^, D), 

D = qq-\^ D » = W ~ + ¥ € ^> (3 ' 15) 

D = 3 + ( 3 - 16 ) 
89 * 

which anti-commute with the supercharges, {Q a ,Db} = 0, and satisfy the algebra, 

{D, D^} = -id+, {D,D^} = ie^d~, {others} = 0. (3.17) 

A chiral superfield $(a;, 6a, 9 a ) and an anti-chiral superfield 9a, 0a) on the lattice 
can be defined in terms of the supercovariant derivatives (13.1 5ft - ( l3~T6l) . 

[SD,${x,9 A ,0 A )]* = [iD,${x,6 A ,6 A )}* = 0, (3.18) 
[^D^{x,6 a ,6a)}* = (ji : no sum), (3.19) 

here and in the following we shall represent the arguments of the superfields by 9 A = 9^ = 
(6\, 62) and 9 a = (9, 9). In solving these conditions, it is convenient to introduce the following 
difference expansion operator, 

where the exponential with the star is defined by 

ef = 1 + X + -X*X + -X*X*X + --- . (3.21) 
2 3! 

The supercovariant derivatives (I3.15ti3.16l) can be expressed in terms of U as, 

= u- l *i^-*u, Id = u- l *i^L*u, (3.22) 

89 89 



^ = U*^^* 11 " 1 no sum), (3.23) 

OK),, 
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with which one may rewrite the chiral and anti-chiral conditions (I3.18l) - (l3.19p as 

[t— ,&(x,e A ,d A )]. = [i-~,&(x,e A ,e A )\. = o, (3.24) 

ott 09 
[^^-^\x,6 A ,e A )l = (ji : no sum), (3.25) 

where 9 A , 9 A ) and $ (x, 9 A , 9 A ) are defined by 

&(x,9 A ,9 A ) = U*<5>(x,9 A ,9 A )*U~\ W(x,9 A ,9 A ) = U' 1 * 9 A ,9 A ) * U. (3.26) 

One can easily solve the conditions (13.241) and (13.251) and get the following expressions for 
$' and $ , 

&(x,e A ,e A ) = &(x,e A ) = tw + e^w + e&fc), (3.27) 

t(x,e A ,9 A ) = t(x,9 A ) = tp(x)+9 X (x) + 9x(x) + 99$(x), (3.28) 

where (0, 0, ip, (p) and (ipi,ip2, Xi X) are denoting the bosonic and fermionic component fields, 
respectively. Then the expressions of the chiral and anti-chiral superfields ($, $) can be found 

as, 

<z>(x,9 A ,9 A ) = ir 1 *^ , (x,e A )*u 



<t>{x) + e^x) + 9 1 9 2 4>{x) + -00 M A M 0(x) - -e^99^A u <p{x) 

+ i e flu 99 1 9 2 A^ u (x) - Uetf^Mx) + -999 ,9^^ ^{x), (3.29) 



®(x,9 A ,9 A ) = U*^'(x ) 6 A )*U- 1 



<p(x) + 9 x (x) + 9~ x (x) + 9~9y(x) - % -99^ip{x) + % -e^99^A vV (x) 



+ t -e^999 fl A uX (x) + % -9~99^x{x) + ^M&A^A^x), (3.30) 

where A M is denoting a symmetric difference operator, A fM f(x) = ^(f(x+n fl )—f(x — n ft )). In 
deriving these expressions we made use of a star product version of the exponential formula, 

e**Y*e; x = Y+[X,Y]* + ^[X,[X,Y].] m + --- , (3.31) 

which obeys from the associativity of the star product. 

The J\f — D — 2 twisted SUSY transformation laws for the component fields follow from 
the definition ( 12.211) . As for the chiral superfield <&(x, 9 A , 6 A ), we have 

6 u $(x,9b,9 b ) = U- 1 *5 u &(x,9 B )*U 

= [£,aQ A , U^ 1 * $'(x, 9 B ) * [/]*, (3.32) 

from which it obeys 

5 u &(x,6 B ) = [UQ A ,&(x,9 B )} # , (3.33) 
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5? 


Six) 





fx) 





ib,, (x) 


— i£ A,, Six) 




ie,, n £ A n Six) 

^ppS — pV^ V / 


Six) 


zfe^A (a;) 

S pO p T fJ \ J 





— i£A n ib n ix) 

S p T P\ J 


rA* t V 




o 


S> A V / 







-i£ p A p ip(x) 








ie pa £, p A a ip(x) 





(p{x) 





i£pA P x(x) + i£ p e pa A a x(x) 






Table 1: N = D = 2 twisted SUSY transformation laws on the lattice for the component fields 

0,Vv><A) and (<p,x,x,<p)- 



where £,aQ'a = U * £,aQa * U 1 are given by 



d , d 

£Q' = f ^ + if0 M dJ, ^Q'p = ^P-Qf (/i : no sum), (3.34) 



lQ' = A-iie^d;. (3.35) 

06 

In a similar way for the anti-chiral superfield <3>(x, 6a, 6 a), we have 

S u W(x,9 A ) = [UQa,^'(x,9 a )U (3.36) 
where £,aQ"a = ^ X * ^aQa * U are given by 

ZQ" = Zffl, ^ = ^gf + i^0d+-i^9e fll/ d- (ji: no sum), (3.37) 

lQ" = A- (3.38) 
86 

The M = D = 2 twisted SUSY transformation laws for the component fields are summarized 
in Table [TJ One can verify that the component- wise SUSY variations form the following off- 
shell closed algebra, 

&>%,W = +i£€pAp(p A (p: no sum), (3.39) 
[S^6^ p )(p A = -iiC p e pa A a cp A (p: no sum), (3.40) 
[others\if A = 0, (3.41) 

where ip A denotes any of the component field (0, ipp, S, ip, x, X, 0)- 
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3.2 Twisted SUSY invariant action at tree level 



In terms of the chiral and anti-chiral superfields <&(x, 9 a) and 9a), one can generally 
construct a M = D = 2 twisted SUSY invariant action on the lattice as 

s = Y.f di ° K *(®> $ ) + Z / ^ + Z / ^ ^® (3 - 42) 

where the Grassmann measures are defined by g? 4 # = d9d9d9id92, d 2 9 = d92d9\ and <i 2 # = 
cf0d0, and each Grassmann integration is defined by J cf# A 1 = 0,/ d9A9s = <^ab- The 
symbol If* is denoting any functions of the chiral and anti-chiral superfields in terms of the 
star product, while the symbols and F* are respectively denoting any functions of the 
chiral superfields and anti-chiral superfields again in terms of the star product. Since we 
take the symmetric choice of a a (13.121) . the summation over x should cover (int, int) sites as 
well as (half-int, half-int) sites, 

= + (mi, m 2 : integers). (3.43) 

x (mi, ma) ( mi + | >m2 +i) 

The necessity of summing over both (int, int) sites and (half-int, half- int) sites stems from 
the fact that in the symmetric choice of a a the lattice SUSY transformation can essentially 
be regarded as a mapping from (int, int) to (half-int, half-int) or vice versa, due to the non- 
commutativity associated with 9a and £4. The simplest action with mass and interaction 
terms can be given by a TV = D = 2 Dirac-Kahler twisted analog of the Wess-Zumino model 

HB1, 

YYl q <yji n 

K* = $ * $, F* = * $ + * $ * $, — * $ + * $ * $. (3.44) 

The total action Swz is accordingly given by, 

SwZ — Skin + S mass + S'jnt, (3.45) 

where the kinetic terms S^m the mass terms S mass and the interaction terms Si nt are ex- 
pressed in terms of the component fields as, 



Sun = Yl Jd , e^(x,e A ,e A )*4>(x,e A ,e A ) 



E 



<f(x)A fl A fl (f)(x) + ix(x)Au,ipu(x) + ie^x(,x)A n ip u (x) + <p(x)4>(x) 



(3.46) 



Sn 



"9" 



d 2 fl $(x, A , A ) * $(x, A) A ) + / d 2 $(x, A , A ) * $(x, A , A 



m 
~2 



E 



(x)(j)(x) + 4>(x)(f)(x) - e^ipJxj^x 



+(p(x)<p(x) + <p(x)<p(x) - x(aOxfc) + x(s)x(x) 



(3.47) 
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s 



int 



-Y 

3! ^ 



d 2 9 $(x, 9 A , 9 A ) * $(x, 9 A , 9 A ) * #4, ^ 



+ / $(x, X , A ) * $(x, 9 A , 9 A ) * A , 9 A ) 



-Y 

3! ^ 



[x)ip u (x + a u )4i u (x - a M ) - e M „W( x + a-v)i>v{x - a u )(f)(x) 



-e^ipfj,(x - a v )4>{x)ijj v (x + a M ) + <fi(x)(fi(x)(p(x) + <f)(x)4)(x)<j)(x) + 4>(x)4>(x)4>(x) 
-(p(x)x(x + a)x(x - a) + (f(x)x(x + a)x(x - a) - x(x + a)x(x - a)<p(x) 
+X(x + a)x(x - a)(p(x) - x(x - a)<p(x)x(x + a) + x(x - a)(p(x)x(x + a) 



+(p(x)(p(x)(p(x) + (p(x)0(x)ip(x) + 0(x)<p(x)ip(x) 



(3.48) 



where we used the symmetric nature of the shift parameters ( 13.121) . a\+a2 = 0, a+a = 0. One 
should notice that the kinetic terms (I3.46[) and mass terms (13.471) consist of the component 
fields living only on (int, int) sites or (half-int, half-int) sites, which implies that the field 
propagations are restricted only from (int, int) to (int, int) or from (half-int, half-int) to 
(half-int, half-int) sites. On the other hand, the interaction terms (13.481) involve both (int, 
int) and (half-int, half-int) sites. In particular, the Yukawa coupling terms consist of the 
bosonic field at an (int, int) site and the fermionic fields at (half-int, half-int) sites or vice 
versa. 

The twisted J\f — D — 2 SUSY invariant nature of the total action is manifest from 
their superfield expressions, since the star products of the superfields satisfy the Leibniz rule 
(I2.27P which is also valid for a product of three superfields, 



5a($i * $ 2 * $2 



6a$1 * $ 2 * $3 + $1 * 5 A &2 * $3 + $ x * $ 2 * 5 A $ 3 . (3.49) 



It is important to note that the above star product nature is implicitly inherited also in the 
component-wise expressions. Actually, in order to be consistent with the star product nature 
of the superfields which preserves the Leibniz rule, one needs to recognize that the SUSY 
variation of each component field should be associated with the star product, for example, 



5 A (<fr(x)(j)(x)4>(x)) 



5a4>{x) * 4>(x)(p(x) + (f>(x) * 5a4>(x) * (p(x) + (p(x)4>(x) * 5a4>{.x). 

(3.50) 



Since the component-wise SUSY variations are defined together with the Grassmann pa- 
rameter £/i's, one needs to take care about the ordering of the component fields in showing 
the SUSY invariance, even though each component field is defined as an ordinary bosonic or 
fermionic function. For instance, 



6 A (<j>(x)4>(x)) ^ 6 A (<i>(x)<i>(x)), 



(3-51) 



even though we have (j)(x)<j>(x) = <f>(x)(f>(x). These properties of the lattice SUSY transfor- 
mation laws essentially necessitates the notion of "proper" ordering in the component-wise 
expressions. For this reason, we respected the orderings of the component fields in the 
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expressions (I3.46I) - (I3.48I) . It is a straightforward task to verify that these proper ordered 
component-wise expressions for the action fl3.46p - fl3.48l) are manifestly invariant under all 
the N = D = 2 twisted SUSY transformations on the lattice. 

We note here that the notion of "proper" ordering in the component-wise expression has 
been already argued in Ref. [H E] in replying to the critiques claimed in Ref [5], and now 
we are revisiting this issue again with the star product formulation at hand. It is obvious 
that the necessity of "proper" ordering stems from the intrinsic non-commutative nature 
of the star product defined in (12.81) and is ultimately originated from the nature of the 
Leibniz rule for the difference operations. Now the question is how this kind of "proper" 
ordering can be survived or respected in the quantum treatment of the formulation. This 
is actually the main subject of this paper. We anticipate the result and claim here that 
respecting the "proper" ordering essentially corresponds to picking up the planar diagrams 
in the perturbative expansions. Developing the superfield formulation, we will explicitly 
see in Sec. 5 that in the planar diagrams all the phase factors which stem from the non- 
commutativity strictly protect the SUSY structure on the lattice and give the exact SUSY 
result. In other words, the above notion of "proper" ordering can be preserved and respected 
in the planar limit or in the 't Hooft large-iV limit. Thus, the exact SUSY w.r.t. all the 
supercharges at non-zero lattice spacing is achieved in this limit. 

We also note some important properties of superfield integrations under the x-summation 
(I3.43p . For the kinetic term or the D-term, we have, 

Y [d 4 9$*$ = Y [d 4 9$$ = Y I rf4 # $ *> ( 3 - 52 ) 

X XX 

which stems from the vanishing sum of the shift parameters (13.111) . For the F-terms in 
particular for the mass terms, we have, 

Y I d 2 9 $! * $ 2 = I d20 $ i $ 2 = Y / d * d $2$1 ' ( 3 ' 53 ) 

X XX 

Y I d 2 d$ l *§ 2 = Y [ d 2 9 $i$2 = Yl I d 2 0$ 2 $i, (3.54) 

X XX 

thanks to the symmetric choice of the shift parameter (13.121) . Namely, the kinetic and 
mass terms can be written without the star product under the x-summation ( 13.431) . This 
is a preferable feature for the superfield formulation, particularly in deriving the superfield 
propagators as we will see in the next subsection. For the interaction terms, we have the 
cyclic permutation property for the (anti-)chiral superfields, 

1 * $ 2 * $ 3 = Yj J d 2 9 $ 2 * $3 * $! = Y2 J d 2 9 $3 * $! * $ 2 , (3.55) 
Y / d 20 ^1 = Y d 2 9 $2 *®3*®i = Y / rf2 ^3 *^2, (3.56) 

X XX 

which is again thanks to the symmetric choice (I3.12p . Actually, one can show that the 
equalities in (I3.53P and (13.551) hold only if one takes ai + a 2 = 0. Likewise, the equalities 
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4- 



--0- 



--Q- 



••0- 



• e (mi,m 2 ) 
■ e (m\ + \,m% + \) 
o g (mi + j, 7712' 



-•— □ e (mi,m 2 + \) 



Fig. 2: Configuration of the superfields in the action subject to the x-summation (|3,57|) 



in (13.541) and (I3.56j) hold only if one takes a + a = 0. This technical advantage is the main 
reason why we take the symmetric choice (13.121) . 

Although the summation defined in (I3.43P is already sufficient for ensuring the SUSY 
invariant nature of the total action, in the following for the technical convenience in the 
momentum analysis, we will conventionally take the x-summation to cover (int, half-int) 
and (half-int, int) sites as well, 

5^ = 5Z + + + ( m i' m 2 : integers). (3.57) 

x (mi,m 2 ) ( mi +i )m2 +i) (mi+|,»na) (roi,m 2 +§) 

Namely, we consider a finer lattice which can be measured by half of the unit length. Ac- 
cordingly, the field propagations are now also allowed from (int, half-int) to (int, half-int) 
and (half-int, int) to (half-int, int) sites. Notice that the lattice action defined by the x- 
summation over (int, int) and (half-int,half-int) is totally independent from the one given by 
(half-int, int) and (int, half-int) in the sense that they are not correlated each other neither 
dynamically nor supersymmetrically. The multiplicity of the superfields introduced so far is 
summarized as follows. The SUSY invariance essentially requires the (super)fields defined 
on both (int, int) and (half-int, half-int) sites as in (13.431) . This causes the field multiplicity 
of factor 2. Furthermore, as explained above, we conventionally introduced a trivial mul- 
tiplicity of factor 2 which stems from the summation (I3.57p . The entire multiplicity of the 
superfields originated from the x-summation is then 

2(SUSY) x 2 (conventional) = 4. (3.58) 

We observe that the lattice action subject to the x-summation given in (13.571) describes four 
copies of the D = M = 2 twisted Wess-Zumino model in the naive continuum limit. We will 
come back this point when dealing with the momentum space analysis of the superfield prop- 
agators in the next section. The lattice superfield configuration subject to the x-summation 
(13371) is depicted in Fig. H 
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3.3 Superfield propagators 



Having defined the M — D — 2 twisted SUSY invariant action on the lattice, we shall 
proceed to set up the technical ingredients requisite for studying the possible quantum 
corrections. We shall take a full advantage of the lattice superfield method and intend 
to calculate the radiative corrections entirely in the superfield method. In what follows we 
derive the superfield propagators on the lattice. The notion of superfield propagators was 
first introduced in [131 [T5]. The derivation in this section basically follows the same procedure 
as in the continuum spacetime [151 ELY] or in the Moyal non-commutative spacetime [19] . We 
start from the free part of the tree level vertex functional, which we shall denote as r . It is 
given by the free part of the lattice action (I3.45p . Using the relations (I3.26p . we express it 
in terms of $' and $' in fl3T27]) - (l3T28|) . 



r [$,$] 



E 

X 

E 



d*eW(u- 2 *&*u 2 ) + ^J2 



J d 2 9 + 


f d 2 d$%' 


J d 2 9 + 


[ d 2 ett 



(3.59) 
(3.60) 



where the U is denoting the difference expansion operator defined in (I3.20p . and the Grass- 
mann integration measures are again defined by d A 8 = d8d8d8id$2, d 2 9 = d02d8i and 
d 2 9 = d6d9. The corresponding connected Green's functional, which we denote as Wq, 
is defined by the Legendre transformation, 



w [j, J] = r [$, $] + J d2 ° J,$/ + J d20 



(3.61) 



where J' and J denote the chiral and anti-chiral source superfields which obey the same 
conditions as $' (13. 24ft and $ (I3.25j) . respectively, 



Q / _ 



[i®j>( x ,e A )l = o, 

ad 



0, fa 



no sum 



(3.62) 
(3.63) 



The Legendre transformation (13.611) gives rise to the relations between (J', J) and ($', $ 

^r T/ 5T 



5$' 



-J', 



5$ 



-1\ 



(3.64) 
(3.65) 



The superfield propagators for $' and $ are defined in terms of the functional derivatives 
of the connected Green's functional w.r.t. the source superfields (J', J), 



< $ (1W(2) > 



5 2 



5 J (1)5 J' (2) 



W [J,J] 



57(1 



-mm, 



(3.67) 
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where the arguments 1 and 2 refer to each superspace point. In the second equalities, we 
used the relations (13.651) . 

In order to calculate the superfield propagators, we first need to write down the l.h.s.'s 
of the relations (I3.64p explicitly, and then solve them w.r.t. $' and $ . By functionally 
differentiating the expressions ( 13.591) and (I3.60f) by $' and $ , respectively, we actually have, 



d 2 9 {1) U 2 (l) * $'(1) * U~ 2 {1) + m<S>'(l) 



J d 2 9 w u- 2 (i) * $'(1) * U 2 (l) + m&(l) 



-J'(l) 



-J (I) 



(3.68) 
(3.69) 



In deriving them, we made use of the functional differential relations for ($', $ 
<5$'(2) 



5 2 xw , 2) 5 2 (ev-ew), 



5$ (1) 



(3.70) 



where 5 2 w (2 ) denotes the two dimensional Kronecker delta, while the delta functions for 



the Grassmann coordinates 5 2 (9^ — 9^) and 5 2 {9^ — 9^) are defined by 



^^(1) _ (2)) = Y[(e%-e™) = (ep-erwF-e™), (3.71) 

A 

S\9 [1) -9 (2) ) = H(9^-9f) = (0U-eW)(p)-p)). (3.72) 



Operating U 2 * from the left and *U 2 from the right on (I3.69f) and taking the d 2 9 integra- 
tion would give 

Km + m J^ *.<!,. tftD.tr-a) = u 2 {i)*7(i)*u- 2 {i), ,3.73) 

and the similar expression for (I3.68P by operating U~ 2 * from the left and * U 2 from the 
right. In (13.731) we used the following relation which can be verified explicitly, 



d 2 9 U 2 * 



d 2 9 U~ 2 *&(x } 9 A )*U< 



*U- 2 = Al<l>'(x,9 A ). 



By substituting (I3.68j) into (I3.73P and solving it w.r.t. $' we then obtain 



$'(1) 



-1 



A 2 — m 2 
t 1 



d 2 9 (1) U 2 (l)*j\l)*U- 2 (l) + 



A 2 — m 2 
t 1 



(3.74) 



J'(l). (3.75) 



Likewise we obtain for $ 



$(1) 



A? ( — m 2 



d 2 Q {l) f/ - 2 (l) * J'(l) * ^(1) + 



m 



A 2 — m 2 



J(l). (3.76) 
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By functionally differentiating the relations fl3.75p - fl3.76l) w.r.t. the source superfields J' and 
J and using the definitions fl3.66p - fl3.67l) . we finally obtain the superfield propagators which 
are given by 



< $VMV(- (2) ,ei 2) ) > = ^^A, x ^\e {1) ~ ^ (2) ) 



2n {Ait) 2 sin p^ + m 2 



A 2 — m 



_ 2 x\ 

" 2T jP- e ^-^ 3 >) ~ m J 2 (e {1) -6 { \ (3.78) 
sin p M + m 2 





(4vr) 2 




1 


^" 


2 e 

- m 2 


/•27T 


d 2 p 


J-27T 


(4tt) 2 




1 




2 e 

- m 2 


/•27T 


d 2 p 


i-27T 


(4vr) 2 



c (l).a.(2) 



sin p M + m 2 



sin p M + m 2 

where the symbols _E (12 ) and are defined by E$ j) = 9^9^ +e^9^9 ( J ] for z,j = 1 or 2. 
In each second line of these expressions, we inserted the momentum space representation of 

2 

c (l) X (2)1 



, n(21) ■ 



r> 2 



"2tt >2„ 

A 2 - / _^_P_Jp{xW-xW) (r>M\ 



Notice that the the range of the Brillouin zone is 47r for each direction since we introduced the 
superfields located on the half-integer sites for each spacetime direction. Within that range 
of the Brillouin zone, each propagator fl3.77l) - fl3.80l) has superficially 16 spectrum doubling, 
namely, 4 species for each direction. 

It is important to recognize that the entire doubling is originated from the two reasons. 
The first one is due to the x-summation over the half-integer sites f)3.57p . This can be 
regarded as a trivial multiplicity as far as the superfield propagators are concerned, since 
the field propagations are only restricted from (int,int) to (int,int), (int,half-int) to (int,half- 
int), (half-int,int) to (half-int,int) and (half-int,half-int) to (half-int,half-int) sites. In the 
configuration space, this is clear from the fact that each term in the kinetic action (I3.46P 
consists of (int,int), (half-int,half-int), (int,half-int) or (half-int,int) sites. In the momentum 
space this feature is reflected to the 2n periodicity of the integrand in each propagator except 
for the exponential factor. In fact, starting from the above expressions of the propagators 
fl3.77p - fl3.80p . one easily sees that there are no correlations between the integer and the 
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half-integer sites for each direction. For instance, the propagator (13.77P can be re-expressed 
as 

r e *«ffl-»M) ~ m 5 2 (0U - gOO), for a;( 1 )-x( 2 ) = (m 1 ,m 2 ) 
i-TT (2vr) 2 sin 2 p M + m 2 (3.82) 

0, otherwise 

where mi and m 2 denote any integer values. Note that the range of the Brillouin zone for 
the non-vanishing propagator is now taken as 2ir since the integrand is actually 2ir periodic. 
In deriving this expression, we divided the momentum regions into two of 2tc and made use 
of the 2ir (anti-)periodicity of the exponential factor e ip ^ (1) '~ x{2) ) . The other propagators 
fl3.79p - fl3.80p also have the same feature. This observation implies that the factor 4 in the 
entire multiplicity due to the x-summation over half-integer sites are just corresponding to 
the four copies of the physically equivalent momentum regions as long as the propagators 
are concerned. 

In contrast, the rest of the spectrum doubling which still remains in the propagators has 
a non-trivial physical meaning via the D = M = 2 Dirac-Kahler mechanism. This stems 
from the fact that the chiral and anti-chiral superfields ($', $ ) contain the components of 
the Dirac-Kahler twisted fermions (XjV'/oX) which turn into the corresponding staggered 
fermion components on the lattice. In this procedure, the extended M = 2 SUSY index can 
be regarded as Nf = 2 flavor (taste) index via the Dirac-Kahler expansion of the component 
fields on the lattice, 

foi(a;) = (lx{x) +'Y^ li (x + ri li ) + r y 5 x(x + n 1 +na))ai, (3-83) 

where the suffices a = 1, 2 and i — 1,2 are representing the spinor and the internal indices 
of N = 2 twisted SUSY. This procedure has been already pointed out in [1] and it resolves 
the multiplicity of 4 into 2 (spinor) x 2 (A/" = 2). In order to see how this mechanism is 
incorporated in the present superfield formulation, it is instructive to derive the propagators 
for the component fields in terms of the £ a ; basis which is given by 

< U(x {1 %p(x {2) ) > = < (lx(a: (1) ) + 7^(^ (1) + n M ) + 7sX(^ (1) + ™i + n 2 )) ai 

x (lx(a: (2) ) + l^M 2) + *V) - 75X(^ (2) + n x + n 2 )) jf3 >,(3.84) 

where £ is defined by the transpose of £, £ = £ T , and the 7 M matrices are given in (13.41) . 
Since £ ai is made up with the component fields located in every corner of the square with 
one integer lattice spacing each side f)3.83p . the locational separation \x^' — x^] in terms 
of £ a j should be taken as even number for each direction. This implies that the physical 
range of the Brillouin zone for £ Qi should be restricted as ir instead of 2ir for each direction. 
Each correlation function for the component fields can be obtained by differentiating the 
corresponding superfield propagator w.r.t. #'s, for instance, 

< x^U^ +»„) > = Jiyi < *(xVo A )#{xV + nil ,e A ) > Cr A 7 even 

(27r) 2 sin 2 pp + m 2 
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Making the use of the completeness relation for the 7^ matrices, 



(l)ij(l)fcj + (ipdki + {lh)ij{l5>ki = 28i k 8ji, (3.86) 
one obtains the following expression of the propagator for the basis, 



_2L (27r) 2 sin p p + m 2 



x 

A' 



(l^apSij sinpp - 2^(75)^(757^)^- sin p M - M a/?) y 



,(3.87) 



where the mass matrix M a p^j is given by 

M aPtij = 2m [<5 a/3 (7 5 )ij cospi cosp 2 + (75)^^ sinpx sinp 2 

-i{li)<xp{l2)n cos Pi sin P2 + i{l2)ap{li)ij sinpi cosp 2 ] • (3.88) 

The integrand in (13.871) can essentially be regarded as the momentum space expression of the 
two dimensional Nf = 2 staggered fermion propagator in terms of spin-flavor(taste) basis. 
Notice that because the range of the Brillouin zone is 7r for each direction, we do not have any 
spectrum doubling in this fermion basis. Also note that thanks to the manifest superfield 
formulation, this mechanism also provides a physical interpretation for the spectrum doubling 
in the bosonic sector. To summarize, the entire multiplicity which is 16 can be decomposed 
into, 

16 = 4(Dirac-Kahler) x 4(x-summation) (3.89) 
= 2(spinor) x 2(W = 2) x 2(SUSY) x 2 (conventional), (3.90) 

where, as we have seen, the factor 4 (x-summation) can be regarded as the trivial multiplicity 
as far as the propagators are concerned, while for the factor 4 (Dirac-Kahler) the Dirac- 
Kahler mechanism plays a physically substantial role. 

Having known the origins of the spectrum doubling, we divide the entire Brillouin zone, 
— 27r < < 2n (/i = 1,2), into 4 x 4 = 16 physical momentum regions by following the 
similar manner as in the staggered fermion momentum space analysis [20] (see Fig. [3]), 

7T 7T 

= + 2(tt k )^ + (-Kl)^ with - - < < -, (3.91) 
where the iix and til respectively with four indices, K,L = 00, 10, 01, 11, are defined by 

ttoo = (0,0), 7Ti = (vr, 0), 7T i = (0, 7r), 7i u = (tt,7t). (3.92) 

In (13.911) , the 2ttk spans the entire momentum space into four copies of physically equivalent 
momentum regions which corresponds to the multiplicity of 4(x-summation), while the itl 
further spans each regions into four smaller regions which corresponds to the multiplicity of 
4 (Dirac-Kahler). The Fourier transformations of the superfields are accordingly given by 
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2tt 



7T 



Fig. 3: 16 physical integration regions in momentum space (solid lines) 



JL2-e**&(p,e A ) = / ^r 2 e l{k+27TK+7TL)x ^K L (k,e A ), (3.93) 

^^'(p,^) = £ / ^e^^^C*,^), (3.94) 

where the superfields in the momentum space with indices K and L are defined by 

& KL (k,9 A ) = &(k + 2Tr K + 7r L ,e A ), t KL (k,9 A ) = t(k + 2n K + 7r L ,e A ). (3.95) 
In terms of the superfields (&' KL , &'kl)-> the propagators are given by 



sin kn + m z 



X -J.7 -/ & m -<?\ (3.97) 
sin + m z 



< ^(i)L(i)(* (1) ,tf )^ao L oo(fc (2) ,fl? ) ) > = (47r) 2 ^ (1))K(2) ^ (1)iL(2) 5 2 (^ 1 ) + fc< 2 >) 

e -4 la) sin^Trx+TTi)^ 1 ) 
X ■ 2,(1)^ 2 » ( 3 - 9§ ) 

sin k/j + m z 

< ^ (1)L(1) (A: (1) ,^ 1) )^( 2 ) L(2 )(A;( 2 ),^ 2) ) > = (^) 2 5 KW , K ^5 LW ,L^ 2 (kW + A;( 2 )) 

+£<? 1) sm(fc+2 7 r if +7r L )« 
X ■ 2, (IK 2 • 

sin kit, + m A 
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Since the momentum range for k is defined from — | to |, these propagators (I3.96p - (l3.99p 
do not have any spectrum doubling in each physical momentum division. Notice that the 
propagators are diagonal w.r.t. the indices K and L. namely, each physical momentum 
division is not mixed up each other while propagating. Also note that as far as the internal 
lines are concerned we are still practically allowed to use the notation p and to integrate 
over the entire momentum region from — 2tt to 2n, keeping in mind that all the 16 physical 
momenta are actually contributing to each internal line. This is because the k, ttk and til 
always come up together and we eventually sum up all the momentum divisions ttk and tyl 
for each internal momentum k. 



3.4 Three-point vertex functions at tree level 

After deriving the superfield propagators, we then turn to consider the superfield three- 
point vertex functions. In contrast to the kinetic and mass terms which can be expressed 
without using the star product (I3.52l) - (l3.54p . the interaction terms cannot be free from the 
non-commutativity stemming from the star product. Namely, the physical implications 
originated from the non-commutativity in this model are essentially inherited in the three- 
point vertex functions. In what follows we derive the superfield three-point vertex functions 
at tree level which we denote and . We start from the tree level vertex functional 

for the interaction terms expressed in terms of and $ in (I3.27l) - (l3.28p . It is given by, 



g_ 

ml gj 



-T 



d 2 e $' + / d 2 e $ * $ * $ 



(3.100) 



In order to properly take into account the effect of the star product, it is the most appropriate 
to Fourier transform the superfields not only to the ordinary momentum space p M but also to 
the momentum counterpart of the Grassmann coordinates 9 a and 6 a which we shall denote 
as Xa = (Xi, X2) and Xa = (x, x), 

&(x,6 A ) = JdpJd 2 xe^e^&(p,XA), (3.101) 
t(x,0 A ) = Jdpjd 2 xe^e^W(p,x A ): (3-102) 

where we used the abbreviations for notational simplicity, J dp = j 2 ^^ t^j-, / d 2 x 

j dx2dxi and / d 2 x = J d%dx while Xa&a = xA + X2O2 and Xa^a = X& + X.8- It is 
understood that the momentum p^ should be decomposed into the physical momenta k^ via 
the relation (I3.9ip whenever necessary. In terms of the expansions (13.1011) and (I3.102p . we 
have for the three-point interaction term in the $' sector, 

j d 2 e $'(x, e A ) * &(x, e A ) * &(x, e A ) 

= \ dH f f[dp {t) d 2 X il) b^'^a * e ^ 2)x e^ * e^*e& 9 * 

x i=l 

x$'(p (1) , X^Wip^lxfWip^, xf), (3.103) 
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and the similar expression for the $ sector. The exponential factors in the r.h.s. of (13. 1Q3[) 
can be evaluated by using the definition of the star product (12. 8|) . 



= e^ XA +Xa +Xa ^ A +p +p ^, (3.104) 

where in the r.h.s. we introduced the notations for the "dressed" xa's defined by, for instance, 
x (3) P u> +P ( 2 > s x (3) e ia A (pV +P W) = x fj(«AU P V+ P V)^ {A . no sum)) etc .. (3.105) 

The is denoting = (01,02) for the $' sector. It is clear that the phase factors in the 
"dressed" form of the Grassmann variables are originated from the non-commutativity in 
the configuration space. After integrating d 2 9 and summing over x and then putting back 
the inverse Fourier transformations, 

$'(p« X W) = [ d 2 6^e- ipil)x( %-*> A &(x^,9%), (3.106) 
x (i) J 

we obtain, 

I d 2 e $'(x, e A ) * $'(x, ^) * $'(x, ^) 

x J 

= EEE / II^V e -^ (1) - (1) +f (2) - C2) +f (3) - C3) ) 52 ^(i) +p(2) +p(3)) 

xW x (2) x (3) J i=l 

jj (£ ,«-p<» -o$){oV*» -ef) ^\e^m^\ef)^\ef), (3.107) 



X 

' A 



where we introduced the notations for the "dressed" #^'s in a similar way as in (13.1051) . 

e^ )p<2> = ^V aAP<2) = Q^e^^, {A : no sum) etc.. (3.108) 

We shall give several technical remarks regarding the expressions (13.1071) and (13.1081) . 
The symbol Sp(...) in (13.1071) is denoting a two-dimensional periodic delta function mod 
An for each direction which obviously stems from the x-summation over integer and half- 
integer sites. Compatibly, the phase factors in the "dressed" 6U's (e.g. ( 13. 108ft ) have also 4ti 
periodicity for each direction. This is clear from the fact that in the symmetric choice each 
a a takes the value of ±| for each direction (I3.12p . The 9 a factors in (13.1071) . 

U(ey-*"-eT) and Ike^-ef) (3.109) 

A A 

are describing the delta functions for the Grassmann variables with the supports, 
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respectively. It is worthwhile to note that each 6 a factor has the following convenient nature 
thanks to the symmetric choice, a\ + a<i = 0, 



IE 



3 (i)-p (3) 



6 



(2)^ 



„(3) 



n^-^ 3 '), (s.iio) 



and in a similar way, we have II a (^a^ > — ^a^) = II a (^a^ ~~ ^a P ' )■ One can use the above 
properties of the 6a factors together with the momentum conservation mod Air wherever they 
are convenient or necessary in the intermediate calculations. 

Let us now derive the three-point vertex functions r$$^ and by functionally dif- 

ferentiating the vertex functional (13.1001) w.r.t. $' and $ , respectively. Using the expression 
(13.1071) . we obtain 



r (tree) ( (1) M) (2) fl (2) (3) fl (3)x 



5 3 r 



r 3 

/ Y[dp {i) e -i(P (1 ^ (1) +P (2) ^ (2) +P (3) ^ (3) ) 52 (pP) + 



p (2) + p (3)> 



X 



n(»r 3 * - s?)^ 



< 3, )+n< 



Of) 



(3.111) 



Likewise for T 



we obtain, 



1 ^ ' P A ' x ' P A ' x J P A / 



5 3 r 



mi 



. 3 

/ e -i(P (1) ^ (1) +P (2) ^ (2) +P (3) ^ (3) ) (pP) + p (2) + 



P 



(3)> 



X 



IE 



,(l)p(3) 



-0 



(3k 



+IE 



j(l)-p( 3 > _( 2 ) _(l)p(2> 7j(3), 
j A ~ U A )\PA ~ U A ) 



(3.112) 



_ -(l)p (3) 
where we introduced the notations for the "dressed" 6U's defined by, for instance, 6 A = 

6^ A e iaAp{i) = 6a e l ( aA ^ pifl ' (A : no sum) with Ha = (a, a). In deriving these expressions, it 
is convenient to use the cyclic properties of the superfields ( I3.55p -( |3~56l) . It is important to 
remark that in the continuum limit, the phase factors in the "dressed" 6U's and 6a§ become 
unity, so that the factors in the square bracket in the vertex functions (13 . 1 1 1 [) and (13. 1 12j) 

are reduced into the ordinary delta functions, 25 2 {6 {1) - 6^)S 2 (6^ - 6^) and 25 2 (d {1) - 
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9 )5 2 (9 —9 ), respectively. On the other hand, on the lattice, the first and the second 
term in each square bracket give rise to the opposite phases each other which essentially 
stems from the non-commutativity between the superfields, $'(1) * $'(2) ^ $'(2) * $'(1), 
etc.. This feature is implying a crucial distinction between the coherent and non-coherent 
multiplications of the phase factors in the perturbative calculations. In the next section, by 
investigating the matrix valued Wess-Zumino model, we will explicitly see that it corresponds 
to the distinction between the planar and non-planar diagrams. 

Before finishing this subsection, we note again that in the above expressions, whenever the 
physical external lines are concerned, the momentum pW should be understood as the phys- 
ical momentum plus corresponding momentum region pW = fcM + 2n K (i) + ir L (e, as defined in 
f)3.9ip . The momentum integration should accordingly be regarded as the physical momen- 
tum integration summed up over all the divisions ttk and 7r L (i), J dp^ = YukW l(») I-i (4^)2 • 
In contrast, as far as the internal lines are concerned, since it K {i) and 7r L (i) always come 
up together and we eventually sum up all the momentum divisions ir K (i) and ir L (i) for each 
internal momentum few, we are practically allowed to use the notation pW and to integrate 
over the entire momentum region from — 2n to 2n, keeping in mind that all the 16 physical 
momenta are actually contributing to each internal line. 



4 J\f = D = 2 Wess-Zumino model with global U(N) 

Since we have derived the propagators and vertex functions, it is now straightforward 
to investigate the possible radiative corrections for the lattice Wess-Zumino model (13.451) . 
However, as we have mentioned in the last subsection, the non-commutativity of the super- 
fields may come up with the crucial distinction between the planar and non-planar diagrams, 
corresponding to the coherent and non-coherent multiplications of the phase factors, respec- 
tively. In order to look at this aspect more carefully, it is appropriate to consider a possible 
matrix extension of the model (I3.45p . In this section, we introduce a N x N matrix version 
of M = D = 2 twisted Wess-Zumino model which possesses U(N) global symmetry. We 
explicitly study the one-loop (in some cases any loops) quantum corrections to the lattice 
action. The model of our interest is given by 



qNxN 

°wz 



qNxN 
kin 



I qlV X iV I C 

' mass ' int 



NxN 



with 



qNxN 
D kin 



O, 



NxN 
mass 



s 



NxN 
int 



/ d 2 9 Tr $ * $ + / d 2 9 Tr $ * $ 

X 

^2 I <,2 9 Tr <I> * <I> * <I> 4- / <P'ti Tr * + * * * 



m 
~2 



9_ 
3! 



(4.1) 

(4.2) 
(4.3) 
(4.4) 



where the matrix (anti-)chiral superfields are defined by <&(x,9,9) = Q a (x,9,9)T a and 
$(x, 9,9) = $ a (x, 9,9)T a with the hermitian generators T a (a = 1, . . . , N 2 ) of U(N). The 
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Fig. 4: Feynman rules for superfield propagators for the U(N) Wess-Zumino model in the matrix 
notation (|4.7p : the momentum should be understood as (k + 2itk + 7I"l)/j if external lines are 
concerned. 



action is obviously invariant under the U(N) global transformation, 

e, e) -> 0, ^jr 1 , e, e) -> 0)fi _1 , with n e c/(jv). (4.5) 

Since all the terms are expressed in terms of the star product, the action is also manifestly 
invariant under the M — D — 2 twisted SUSY transformations, 

8 A $ a = [UQa,**]*, (4-6) 

where the supercharges Qa's are given by (I3.13p - fl3.14p . Deriving the component- wise SUSY 
transformations for the matrix superfields is just in a straightforward manner. Notice that 
the cyclic permutation property of the superfields under the trace is compatible with the 
cyclic nature of the star product listed in f l3.52p -f l3.56p . One can follow exactly the same 
procedures as given in the last section, and derive the corresponding Feynman rules for the 
matrix action (14.11) . If one employs the matrix index notation of the superfields, 

= $ a (T a )>, $/ = $ a (T a )>, (4.7) 
the U(N) matrix extension of the functional derivative relations (I3.70p are expressed as 

Sfl = m^ne {i) -o^i |^ = «^ C2) ^ (i) -^ (2) ). (4.8) 

\ l ) 09 tl (1) 

The superfield propagators and the tree level vertex functions for the J\f — D — 2 Wess- 
Zumino model with global U(N) are summarized in Fig. H] and [5J 

Let us first calculate the one-loop self-energy diagram shown in Fig. [6] which we denote 
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Fig. 5: Feynman rules for the superfield vertices for the U(N) Wess-Zumino model in the matrix 
notation (|4.7|) : the momentum should be understood as (& + 27r K -(.) + ^ L (i))^ if external lines 
are concerned. 



T^ l ° . Up to some overall numerical factors, it is given by 



ioop) 



^ 2 E E/^/ d 2 e®d¥ A) d 2 e (5) d 2 9^ 



+ c c c n< 



3(3) ^ 



X 



<J? B <p 5 i4 TT 
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+ c c c 



7T (2)-g _< 4 ) w7 K2)-(p+9) 77(5) 
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(4). /?T (2)p+ 9 
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B 
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(5)v 
B > 



B 



P +B^ 43) sin(p+g) M -B^ 56) sing M 

* sin (p + qj/j, + m / 



sin 2 + m 2 



(4.9) 



(1— loop) 



where again = 9^9^ + e^O^Qv^ and J dq = jjjjs- After contracting the matrix 

indices, one finds that the one-loop self-energy T 
planar part (Fig. [7]) and the non-planar part (Fig. 



can be divided into two parts: the 



r 



(1— loop) p(l— ioop) -p(l— /oop) 

<J><1> &&(planar) ^^(non— planar)' 



(4.10) 
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Fig. 6: One-loop correction to the vertex function T 





Fig. 7: Planar contribution to T 



(l—loop) 



Fig. 8: Non-planar contribution to T 



(l—loop) 
$5 



which are respectively given by 



loop) 
§&(planar) 



5 n 5 n S?N f dq f d VV£ (4 W 5) d¥ 6 ) 



7j(2) ~q 7i( 4 ) \ /7)( 2 ) - (P+9) _ fl( 5 ) n 
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+ U( e A P+q - of)(e^ - of) H(9% q - ef)(ef p+q - ef 

A B 
e +Bp 43) sin(p+ij) M g-Bp j6) sinij M 
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sin {p + q) n + m 2 sin g M + m 2 



(4.11) 



-p(l-ioop) 

<E><E>(non— planar) 



gj; S fl 9 - I dq I dW5 (4 WV9< 6 > 



n< 



< 6, )(»! 



3(2)9 

1 B 
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ii 



3 +E< 43) sin(p+ 9 ) p 



X 



c(56) ■ 
e -E^ 'smij M 



sin (p + q) u + m 2 sin g M + m 2 



(4.12) 



Note that the planar part has a factor iV which stems from the trace of the matrix indices. 
One can define the 't Hooft large- N limit [21] by rescaling the coupling constant g 2 , namely, 
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by taking iV — > oo with A = g 2 N fixed. In this large- N region, the non-planar contribution 
(I4.12p is suppressed by O(j^) while the planar contribution (14. lip is order of 0(1). 

We shall first evaluate the planar contribution T^ 1 ""^ .. After performing the 6 inte- 

r $<&(planar) r ° 

gration, the first term in the square bracket in (14. lip gives rise to the following argument in 
the exponential, 

[0(2)-ff0(,l)-5 + e ^fif(2)-ff0(l)-«] S in(p + g)^ - ^(2)-(p+ q ) 6 (l)-( P +g) + £ ^gf(2)-(p+«) (l)-(p+«)] gin ^ 

= [^ (2) ^ 1} e-^ + e^^e 4 ^] sin(p + g) M - [fl^We-^^ + e /il i (2) ^ 1) e +i(p+( ^] sing M 
= ^ (2) ^ 1) + ^^ 2) e ) ]sin^ 

= Ef)sin^, (4.13) 

where we used the Leibniz rule conditions (13.91) from the first to the second line. It is 
straightforward to show that the second term in (14. lip also gives rise to the same factor 
after the Leibniz rule conditions are imposed. The planar contribution to the one-loop 
self-energy diagram T^ loopS> is thus given by, 

r (l-ioop) = 9 2 N j j +E w sinpu f 2n d 2 q 1 1 (4U \ 

^(planar) 2 11 12 J _ 2n (4tt) 2 sin 2 (p + q) ^ + m 2 sin 2 q^ + m 2 ' l ' ' 

Since the superspace coordinates are exclusively encoded only in the exponential factor 

(21) 

e E i* smPf * in (14.141) . the M = D = 2 twisted SUSY invariant nature of the kinetic part of the 
action (13.461) is strictly protected at non-zero lattice spacing as long as the planar diagram 
is concerned. In other words, the supercoordinate structure and matrix index structure in 
the superfield propag ator < $'(1)^$ (2)<f > remains unchanged even after the planar one- 
loop contribution is taken into account. The numerical factor in (14.141) which stems from 
the g-integration may be absorbed into a wave function renormalization. These procedures 
are completely analogous to the ones in the continuum theory [TH [16j [17] . Notice that these 
remarkable features are originated from the coherent multiplications of the phase factors 
e iaAP in the planar diagram as shown in (I4.13p . In contrast, as one can imagine, the non- 
planar diagram T^— 1 . 00 ^ , , does not have these favorable features due to the non-coherent 

" <P<p(no?i— planar ) 

multiplications of the phase factors. Actually the first term of (14.121) would give rise to the 
argument of the exponential as, 

10(2)0 9 (l)-a + e ^(2) 9 0(l)- ?] gin(p + g ^ _ [e (2) P+qe (l)- {p+ a) + ^J(2) P+ a e (l)-( P+g)] ^ % 

= sinp M + 2ia ■ (q sin(p + q) n - (p + q) sinp„)(0< 2 W - e^B®) + 0(a 2 ),(4.15) 
while the second term of (I4.12p yields the argument with the opposite phases, 

[0(2)-ff0(l)9 + V 0(2)- 9 0(1) ?] gin(p + q)fi _ [ 0(2)-( P+? )0(l) P+g + ^J(2)-( P+g ) 9 (l) P+q] gin ^ 

= sinp M - 2m • (q sin(p + q)„ - (p + q) sinp /i )(^ 2 ^« - e^O^) + 0(a 2 ).(4.16) 

Although the 0(a) contributions are canceling each other, the non-planar diagram would 
generically spoil the SUSY structure at 0(a 2 ). In the 't Hooft large- iV region, however, 
these SUSY breaking contributions in r^ 1 -^ 00 ^ , , are suppressed by 0(a 2 /N), since the 

° <&<&(non— P lanar) ^ L J K 1 11 
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Fig. 9: One-loop correction to the vertex function 





Fig. 10: Planar contribution to T^ loop ^ 



Fig. 11: Non-planar contribution to T^ loop ^ 



coupling constant g 2 scales as jj. We thus state that at one- loop level, M — D — 2 twisted 

SUSY is strictly protected at non-zero lattice spacing in the self-energy vertex T^ loop ^ only 
in the planar limit, or equivalently, in the 't Hooft large- N limit. 

We next evaluate the one-loop correction to the vertex function (Fig. [9]) which 
given by 



is 



-p(l-Zoop) 
1 <J>$ 



g 2 J dqf d 2 e^d 2 e^d 2 e^d 2 e^ 
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x sit S h 



2 3 "l 4 



-m(5 2 (^ 3 ) -0( 4 >) 
sin 2 (p + g) M + m 2 '"' ,( 



sin g M + m 2 

Following the same procedure as in T^ loop \ we can divide the vertex diagram r$$ loop ^ i 
two parts: the planar diagram (Fig. [TOl and the non-planar diagram (Fig. [TTT) . 



(4.17) 



-,(1— loop) 



■p(l-Zoop) -p(l-Zoop) 
$$(p/anar) <l?<I>(rion— planar) ' 



(4.18) 
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Fig. 12: One-loop correction to the vertex function 



After performing the 6 integrations, one finds that integrand in the planar diagrams is 
proportional to, 



= iW - - <>f) + im 11 - - «!?> 

= 0, (4.19) 

where we used the relation ai + ai = and factored out the overall phases along the 
same procedure as in (13.1 10p . We thus have the vanishing contribution from the planar 
diagram. Note that this planar calculation is completely analogous to the calculation in 
the continuum spacetime. Namely, loop ^ should vanish if the exact SUSY invariance was 
really maintained. In contrast, the integrand in the non-planar diagram is proportional to, 
after the 9 integrations, 

1 ^(non-planar) K [[^A A )[?A ~ d A ) + [[^A A W A ~ 6 A 

A A 

= 8^ 1) ^ 2) 4 1) ^ 2) sin(ai-p) sin(a 2 -p) ~ 0(a 2 ), (4. 

which implies that the non-planar contribution generally spoils the exact SUSY structure 
for the non-zero external momentum p. In the large- iV region, however, this contribution 
is again suppressed by 0(a 2 /N) as in the case of T^ 1 . 00 ^ , .. In a similar way, we can 

° 1 r ^^(non— planar) J1 

show that the planar contribution to the vertex function of the anti-chiral superfield T^ loop ^ 
manifestly vanish with the same mechanism as the above, while the non-planar diagram 
generically spoils the SUSY structure at 0(a 2 /N). 

Let us turn to the one-loop correction to the 3-point vertex function T^q,® (Fig. [12]) 
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Fig. 13: Planar contribution to T^^° p ^ 



Fig. 14: Non-planar contribution to rj^^ 00 ^: 
we also have the diagrams with cyclic permuta- 
tion of (h,ji), (12,32), (h,h)- 



which is given by 
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(4.21) 
(4.22) 



Following the same procedure as before, one can see that the planar contribution r$$$°p Zanar ) 
(Fig. [T3|) . which stems from the coherent multiplications of the phase factors turns out to 
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be vanishing. Actually, after the 9^ l '(i = 4 ~ 9) are integrated, it is proportional to, 

1 ^(planar) x H^A _ ^A "A )\ U A ~ U A ) 

A 

A 

A 

+n(«i i, "-«r)(«i i, "-«f')(«?"-«r) 

A 

= 0, (4.23) 

where from the first to the second line, we made use of the properties of the 9a factors such 
as (13. 110p . while from the second to the third line in (14.231) . we absorbed the phase factors 
and redefined the 9a's, 

9f = 9^~\ 9f = ef^ W+p ^\ of = 9f~ {p{2)+(1 \ (4.24) 

e (iy = e o.* e w> = e w*^« e w = qww^ (425) 

in order to emphasize the manifest nilpotency of 0® 's and 9® 's. On the other hand, the 
non-planar contribution ^^^"non-pianar) (Fig- D21) turns to be proportional to 

^Lptanar) * 8 9? 9? 9^ 9? 9? 9® [sin a, • ( Pl + p 2 ) sin a 2 • ( Pl + p 2 ) 

+ sin ai ■ pi sin a 2 • pi + sin a\ • p 2 sin a 2 • p 2 ] 
~ C(a 2 ), (4.26) 

which is not vanishing for the generic external momenta p^ and p^ . In the large- N re- 
gion where the coupling constant g is proportional to the SUSY breaking non-planar 

contributions to the tree point vertex function rj^ oop ^ are thus suppressed by 0(a 2 /N^). 
In a similar manner, one can show that the planar diagrams contributing to the one-loop 
three-point vertex function r^-^ op manifestly vanish at non-zero lattice spacings due to 

the same mechanism as r$$$°p^ nar ), while the non-planar contributions generically spoils 

the SUSY structure at the order of 0(a 2 /N^) in the large- N region just as in the case of 

■p(l — loop) 

<E><J> < I > (non— planar) ' 

Let us summarize the above results for the superfield one-loop calculations at non-zero 
lattice spacing. We have seen that the mass m and coupling constant g in the twisted 
TV = D = 2 Wess-Zumino model are strictly protected from the radiative corrections in the 
planar diagrams apart from the wave function renormalization which stems from the self- 
energy diagram f)4.14p . On the other hand, the non-planar diagrams in general give rise to the 



33 



3i H 



/<2) 







9 






2 










(l)i( 3 ) 



? (3)l«< 



(l)-2( 3 > /i(2)-Z( 1 )x/ /1 (l)-K 3 ) 
— (7 



-e 



(3)-Z( 2 )< 




SUSY spoiling contributions although they are suppressed at least at the order of 0(a 2 /N) 
in the 't Hooft large- N region. The breakdown of twisted SUSY in the non-planar diagrams 
is originated from the non-coherent multiplication of the phase factors, which stems from 
the intrinsic non- commutative nature of the star product, $i*$2^2*^i- These features 
imply that the exact lattice realization of M — D — 2 twisted SUSY in terms of the present 
non-commutative framework is achieved in the planar limit, or equivalently, in the 't Hooft 
large- N limit. It should be emphasized that, as we have seen in the above calculations, 
the loop calculations in the planar diagrams are completely analogous to the ones in the 
continuum superspace calculations. Therefore, we observe that, by further developing the 
lattice superfield techniques and following the argument of Grisaru, Rocek and Siegel [22] , 
a Dirac-Kahler twisted M = D = 2 analog of the non-renormalization theorem may also be 
proven directly on the lattice in the large- N limit. The result will be given elsewhere. 

Before finishing this section, it is instructive to show within the current superfield tech- 
niques that a certain class of planar diagrams manifestly vanish at any order of perturba- 
tion theory at non-zero lattice spacing. One can actually see that diagrams for (r$$) 
and r$$$ (r$$$) manifestly vanish if the outer edges of the diagrams only contain planar 
vertices and the propagators < > (< <3>$ >). In showing such properties, it is con- 
venient to employ the double-line notation with the conserved momenta fW's defined by, 
pM = /(!) - /(2) )P ( 2 ) = /( 2 ) - /( 3 ) and p (3) = Z (3) - Fig. Q2] depicts the three-point vertex 
function in this momentum notation. Fig. [16] depicts the class of manifestly vanishing 
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Fig. 17: Manifestly vanishing diagram for the three point vertex function r<j,$<j> 



diagrams for the two-point vertex function r$$, which is essentially proportional to, after 
performing the 6 integrations for the propagators, 

(Figna k j d 2 e (1, d 2 e m ■ ■ ■ <w> dV^'dV 2 ' ■ ■ • d 2 9'< m > 



L A 



j(3)Z( 2 >- 



X • • • X 



wr a 



'Wh/n'W (2) /i / (2)«' (1) \//i'(2)«' (3) /i'(3)Z' (2) 

" (7 



+ (the term with the opposite phase) 



B 



0. 



(4.27) 



Likewise the diagrams in Fig. [T7] for r$$$ is proportional to, after integrating the #'s 
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Fig. 18: Examples of the vanishing diagrams with (a) non-planarity and (b) propagators < > 



associated with the propagators, 



:FigH3 « J d 2 e' l >d 2 e^> ■ • • d 2 9'"> dvewro . . . d 2 e'' m > dV' 1 W« ■ ■ ■ d 2 ff" r > 



n< 

• A 



-e 



(3)Z( 2 )< 



X ■ ■ • X 



(n-lWW /) (n)«("- 1 ) n / n (n)V 
— U " " 



,1(2)1 



'(3) 



'(2), 



B 

x . . . x (0j m_1)i " (1) - 6iJ m ) r(m ~ 1) )(6)J m)i " - 



c 



x---x(^- 1) ' (1) -^ )r(r " 1, )(^-^) 



+ (the term with the opposite phase) 



IS - *£) 



0. 



(4.28) 



It is worthwhile to note that the vanishing diagrams Fig. [TS] and Fig. [T7] generally include 
not only planar diagrams but also a certain class of non-planar diagrams whose non-planarity 
reside inside of the diagrams such as Fig. [18] (a). Also the diagrams accommodating the 
propagators < $$ > in their inside such as Fig. 1181(b) also generically belong to the category 
of manifestly vanishing diagrams. 



5 Summary & Discussions 

Having constructed a novel star product formulation of Af = D = 2 twisted SUSY 
invariant action on a lattice, we perturbatively studied the possible quantum corrections 
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at non-zero lattice spacing. We have explicitly calculated the one-loop (in some case any 
loops) radiative corrections for M = D = 2 twisted Wess-Zumino model with global U(N) 
symmetry, and we have seen that the planar diagrams strictly protect the entire M = D = 2 
twisted SUSY structure at non-zero lattice spacing. On the other hand, the non-planar 
contributions generically spoil the lattice SUSY although they are suppressed at least by 
the order of 0(a 2 /N). From these features we claim that the full realization of exact lattice 
SUSY is achieved in the 't Hooft large- iV limit. 

In the series of papers in the collaboration with A. D'Adda, I. Kanamori and N. Kawamoto 
[U [2j [3], we have been proposing a theoretical framework to realize exact SUSY on a lattice. 
Since the lattice action in this framework is given in terms of the superfields, their SUSY 
invariance is formally kept manifest at the classical level, although its quantum behaviour 
has not been clarified until this paper. This paper addresses how the quantum corrections 
would spoil or protect the classical SUSY invariant nature of the lattice action by explicitly 
calculating the radiative corrections to the lattice action perturbatively. As was given in 
the one-loop calculations, (I4.14[) . (I4.19p and (14.231) . the planar diagrams strictly protect the 
SUSY invariance, although the non-planar diagrams in general spoil the lattice SUSY at non- 
zero lattice spacing. This result shows that even though we start from the "manifest" SUSY 
formulation at classical level, the SUSY invariant nature of the action is generally spoiled 
by the non-planar loop conrtibutions. Thus the validity of the manifest lattice SUSY formu- 
lation could be claimed only in the 't Hooft large- N limit at least in the non-commutative 
product formulation introduced in this paper. Nevertheless, it is still important to stress 
that the SUSY protecting nature of the planar diagrams is far from accidental and is actually 
completely analogous to the twisted version of the superfield calculations in the continuum 
spacetime |14|-[18|. This is mainly because the non-trivial phase factors associated with 
each vertex, which stem from the "mild" non-commutativity (12.51) . cancel with each other 
in the planar diagrams. The generally vanishing planar diagrams (I4.27P and (14.281) exhibit 
this aspect explicitly. 

The consequence of this paper may naturally be understood if one reminds the notion 
of "proper" ordering appeared in the following two physically independent contexts. The 
first one is, as we have mentioned in Sec. 3 and also discussed in Ref. [3], that the notion 
of "proper" ordering of (super)fields is essential in the present framework of exact lattice 
SUSY due to the non-commutative nature of the SUSY transformation on the lattice. This 
is originated from the superficial ordering sensitive nature of the difference operations on the 
lattice. The second one is, as is well known, that taking the large-iV limit also essentially 
restricts the ordering of the (super)fields in a certain manner, since interchanging the order 
of the fields generally corresponds to converting planar diagrams to the corresponding non- 
planar diagrams or vice versa. In this paper we have explicitly seen that respecting the 
"proper" ordering in the lattice SUSY context essentially corresponds to picking up the 
planar diagrams in the perturbative expansions. In other words, the somewhat peculiar 
notion of "proper" ordering introduced in the context of lattice SUSY realization now turns 
to have a physically relevant meaning in the large- N limit. 

One might naively think that the large- N reduction [23J would make the lattice SUSY 
problem completely trivial from the beginning since in the large-iV limit the entire lattice 
coordinate may be reduced into a single site, and no chance of the lattice Leibniz rule 
problem itself may occur. However, one should notice that the shift operations or difference 
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operations on the lattice still remain encoded in the shifting matrices even after the reduction, 
for instance, in with V^T^, = Z U ^T U T^, if one takes the twisted Eguchi-Kawai reduction 
|24j . From this point of view, the lattice SUSY formulation in the large- iV limit may naturally 
be converted into a problem of how to supersymmetrize these shift matrices. We should note 
that this is actually the subject which was partially addressed in Ref. [6]. 

It is interesting to consider if the star product formulation introduced in this paper could 
be applied to formulating the supersymmetric gauge theories on the lattice. We observe 
that the lattice gauge covariant formulation introduced in Ref. [21 [3] by means of the link 
supercharges and the link component fields may be extended to a certain type of star gauge 
covariant superfield formulation on the lattice. Also the lattice Chern-Simons formulation 
based on the M = 4 D = 3 twisted SUSY which was recently introduced in Ref. [12] may 
be recasted into this category. It is also worthwhile to further develop the lattice superfield 
framework along the similar manner as given in Ref.|22j and try to accomplish the lattice 
analog of the non-renormalization theorem directly at non-zero lattice spacing. 

In this paper, all the perturbative calculations are performed with the symmetric choice 
of a a (13.121) . We have introduced four copies of the lattice superfields essentially due to this 
symmetric choice. This is because the lattice SUSY transformations involve the mapping of 
the component fields from the integer sites to the half-integer sites, or vice versa, due to the 
non-commutativity associated with 9 a and £4. One may wonder what would happen if one 
takes any other choices. In particular, the asymmetric choice where the shift parameters are 
given by, a = 0, = and a = —n\ — n 2 , would be of the most interest among them. 
As we have presented in Ref. [I], the lattice SUSY invariant action with the asymmetric 
choice can be constructed only on the integer sites, so that one may expect the results 
without introducing any copy degrees of freedom. On the other hand, it is important to 
recongize that the superfield structure, particularly the structure of the vertex functions, 
are strongly related with the rotational symmetry in the configuration space. In particular, 
the relations fl3.53p - fl3.54p and the cyclic permutation relations f)3.55l) - fl3.56p are satisfied 
only if one takes the symmetric choice of a a f l3.12p . In contrast, if we take the asymmetric 
parameter choice, we do not have any cyclic permutational symmetry of the superfields at 
each vertex. In the configuration space, this corresponds to the fact that we do not have any 
non-trivial rotational symmetry if we take a = 0. The lack of the rotational symmetry in the 
case of a = makes the perturbative calculations much more complicated and non-trivial 
even in terms of the superfields. Actually each three point vertex would give rise to six 
independent terms, instead of two in the case of the symmetric choice. Accordingly, their 
N dependences could get less transparent. The perturbative study of the asymmetric choice 
therefore should require more extensive and careful calculations and observations. We will 
keep this subject for future study and the result will be given elsewhere. We should also note 
that addressing the rotational symmetry of this formulation may involve how to define the 
rotational symmetry itself within the context of non-commutative superspace framework. 
We observe that the "twisted" deformation of Lorentz symmetry which has been proposed 
in the context of the Moyal non- commutative spacetime in Ref. [25] may play an important 
role also in the non-commutative lattice SUSY formulation. 

Recently, the Leibniz rule on the lattice was investigated from an axiomatic point of 
view, and its physical relevance with the infinite flavor d.o.f. was pointed out It may be 
worthwhile to look at the non- commutative superspace formulation presented in this paper 
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from that point of view. 
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